Nonlinear viscoelastic dynamics of nano-confined water 
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The viscoelastic dynamics of nano-confined water is studied by means of atomic force microscopy 
(AFM). We observe a nonlinear viscoelastic behavior remarkably similar to that widely observed 
in metastable complex fluids. We show that the origin of the measured nonlinear viscoelasticity in 
nano-confined water is a strain rate dependent relaxation time and slow dynamics. By measuring 
the viscoelastic modulus at different frequencies and strains, we find that the intrinsic relaxation 
time of nano-confined water is in the range 0.1 — 0.0001 s, orders of magnitude longer than that 
of bulk water, and comparable to the dielectric relaxation time measured in supercooled water at 
170 - 210 K. 
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Confined fluids exhibit unique structural, dynamical, 
electrokinetic, and mechanical properties that are differ- 
ent from those of the bulk [J 0, S, i, H, 0, Q- Their 
behavior depends on the degree of confinement, strain 
rate, temperature, fluid molecular structure, and interac- 
tions with boundaries. Surprising effects have been found 
when water is confined in nanogaps. For example, the 
electric field induced freezing of water at room temper- 
ature Q and the extremely high viscosity of water close 
to a mica surface [H, 0, Q . Previous experiments and cal- 
culations have pointed out the key role of the confining 
surfaces 0, 0j- A notable increase in viscosity and de- 
crease in the diffusion constant was measured only when 
water was confined between hydrophilic surfaces. For hy- 
drophobic confinement, the observed increase of viscosity 
was not very pronounced. Intriguingly, a similar behav- 
ior has been observed in confined glassy materials. When 
a glass-forming fluid is cooled down to the glass transi- 
tion temperature, T g , its viscosity grows by many orders 
of magnitude, and the confinement can increase or de- 
crease T g for strong or weak interactions with the walls, 
respectively Q. 

So far, the viscosity measurements for nano-confined 
water have been performed in the linear viscoelastic 
regime. However, as observed in macroscopic rheological 
measurements, the study of the viscoelastic properties as 
a function of shear amplitude and rate is important for 
a better understan ding of the dynamical and structural 
properties of fluids [10( . 

In this letter, we investigate the viscoelastic response 
of nano-confined water, as a function of shear ampli- 
tude and rate, by means of direct high-resolution AFM 
measurements. We observe a nonlinear viscoelastic be- 
havior remarkably similar to that widely observed in 
metastable complex fluids, such as gels and supercooled 
liquids , lllj . We show that the origin of this nonlin- 
ear viscoelasticity in nano-confined water is a strain rate 
dependent relaxation time and slow dynamics. By mea- 
suring the viscoelastic modulus at different frequencies 
and strains, we find that the intrinsic relaxation time, 



To, of nano-confined water is in the range 0.1 — 0.0001 s, 
orders of magnitude longer than that of bulk water, and 
comparable to the dielectric relaxation time measured in 
supercooled water at 170 — 210 K. 

In our AFM experiments [l], [lj| , a nano-size spherical 
silicon tip is brought quasi-statically to the vicinity of a 
flat freshly tape-cleaved hydrophilic mica surface, all im- 
mersed in purified water, while small lateral oscillations 
are applied to the cantilever support [l[ . The normal and 
lateral forces acting on the tip are measured directly and 
simultaneously as a function of the water film thickness, 
i.e., tip-sample distance, d. The zero distance, d — 0, is 
where the normal force diverges. 

The experiments were performed with a Molecular 
Imaging PicoPlus AFM. An instrumental problem in 
quasi-static force measurements is that, during the tip- 
sample approach, the tip snaps into contact with the 
surface at a distance where the gradient of the tip- 
sample forces exceeds the cantilever normal spring con- 
stant, kjsi [3, EH- To overcome this problem, we em- 
ployed relatively stiff cantilevers. We used silicon tips 
with radii R = 40 ± 10 nm and Ultrasharp NSC12/50 
cantilevers with normal and lateral spring constant in 
the range k N = 3 - 4.5 N/m and k L = 50 - 120 N/m, 
respectively. Due to the mechanical stability of our appa- 
ratus, and a judicious choice of the cantilever stiffness [lj, 
we are able to measure the tip-surface distance with sub- 
Angstrom resolution all the way down to the last ad- 
sorbed water layer. The calibration and force detection 
was performed as described in 13, [l6j]- The approach 
velocity was 0.2 nm/s. During the approach, lateral os- 
cillations parallel to the mica surface were applied to the 
cantilever holder by means of a lock-in amplifier. The 
same lock-in amplifier was then used to measure the am- 
plitude of the lateral force, Fl, and the phase difference, 
9, between the applied lateral displacement and the de- 
tected lateral force. The zero of phase was chosen when 
the tip was in hard contact with the mica surface, for 
lateral oscillation amplitudes, Xq, small en oug h to guar- 
antee an elastic contact without slippage [IT] . In order 
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to shear parallel to the mica surface, before each mea- 
surement we tilted the stage that holds the sample until 
the differences in height of the mica surface topography 
across an area of 1 x 1 /im 2 (as obtained from AFM sam- 
ple topography imaging) were smaller than 1 nm. This 
corresponds to an angle < 0.06° between the mica surface 
and the tip shearing. 

All the experiments were performed at 300 K in high 
purity DIUF water from Fisher Chemicals (pH = 6.1). 
The purity of water used in our AFM liquid cell was 
tested after the experiment by gas chromatography - 
mass spectrometry (GC-MS). GC-MS spectra of used, 
and not previously used, water samples were taken by 
70SE spectrometer (VG Instruments). In both cases the 
results showed that any small molecular weight (less than 
700 Da) organic contaminants were present at amounts 
below the instrumental threshold (5 ppm). 

When a viscoelastic material is confined between two 
parallel plates separated by d, with area A, and a sinu- 
soidal strain is applied to one of the plates at the fre- 
quency u>, 7 = 7osin(cj£), the resulting stress between 
the plates can be written as a = <7o sin(cjt + 9). The re- 
lationship between the strain amplitude, 70 = ^f, and 
the stress amplitude, oq = ti is given by the following 
equation: 



(i) 



where G* is the viscoelastic modulus. The viscoelastic 
modulus contains information about the dissipative and 
elastic response of the confined material. In particular, 
G* can be written as a complex sum of the storage mod- 
ulus, G\ and the loss modulus, G", i.e., G* — G' + iG", 
as 
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For a purely elastic solid, a and 7 remain in phase, 9 = 0, 
and so G" = and G' = G* . 

In order to study the viscoelastic behavior of nano- 
confined water we have measured the lateral force and 
the phase when we oscillate laterally the AFM cantilever 
holder. As a first approximation we have assumed that 
the lateral spring constant of our silicon cantilever is 
much larger than the lateral tip-water contact stiffness 
for d < 1 nm fl9j| . As a consequence, we can consider that 
the applied oscillation amplitude to the cantilever holder 
is equal to the shear amplitude of the tip apex. Figure 1 
shows Fl and 9 as a function of d for three different shear 
amplitudes at a fixed shear frequency, ui = 955.3 Hz. For 
tip-sample distances larger than 1 nm we observe that 
the lateral force is equal to zero within the instrumen- 
tal error for any shear amplitude. As soon as d < 1 nm 
we observe that Fl increases with decreasing d, and di- 
verges at d = nm when the tip is in hard contact with 
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FIG. 1: The lateral force and phase as a function of tip- 
sample distance at constant shear frequency, 955.3 Hz, and 
for three different shear amplitudes, (a) Xo = 0.4 nm. (b) 
Xq = 0.66 nm. (c) Xo = 1.32 nm. The phase for d > 1 nm is 
not shown because it largely fluctuates in a random way. 



the mica surface. In a previous study [l|, Fl has been 
used to calculate the viscosity of water (77) by using Eq. [TJ 
and by considering water as purely viscous, that is, by 
making the approximation G* |~ G" ~ r\ ■ u>. This 
approximation is true when 9 = 90°, which, as we show 
later, is the case for large strain rate amplitudes defined 
as 70 = 70 • w. However, the phase measurements pre- 
sented in Fig. 1 show that in general the behavior of 
nano-confined water is viscoelastic, and furthermore, we 
observe that the lateral force does not grow proportion- 
ally with the shear amplitude, nor with the shear fre- 
quency (not shown here). This is an indication that the 
viscoelastic response is not linear, and the viscoelastic 
modulus is shear amplitude dependent, G* = G*(7o). 
Therefore, a detailed study of G* as a function of 70 is 
needed to shed light into this nonlinear behavior. 

By applying Eq. [2] to the data in Fig. 1, we have ex- 
tracted G' and G" as a function of d for different Xo at a 
fixed uj. (The A used for Eq. [2] is the contact area corre- 
sponding to the spherical segment defined by the intersec- 
tion between the spherical tip and a plane at z = d+lXh, 
Ah = 0.25 nm, i.e., a water molecule diameter [lj].) Fig- 
ure 2 shows very clearly that G 1 and G" strongly depend 
on the shear amplitude. G" dominates over G' for large 
shear amplitudes, where the response of nano-confined 
water becomes purely viscous. Also, by decreasing the 
gap size, we observe that the rise of G' and G" takes 
place later (smaller d) for larger shear amplitudes. Fur- 
thermore, for all the investigated shear amplitudes, the 
rise of G" occurs earlier (larger d) than the onset of G' . 
The dramatic drop of both G 1 and G" for d < 0.2 nm 
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FIG. 2: G' and G" as a function of tip-sample distance. The 
shadowed area, d < 0.2 nm, is not discussed in this letter 
because the gap size is smaller than a water molecular di- 
mension. The frequency is 955.3 Hz, and the shear amplitude 
is 0.4 nm for (a') and (a"), 0.66 nm for (b') and (b"), and 
1.32 nm for (c') and (c"). 




FIG. 3: At a constant distance, d = 0.4 nm, G' and G" as a 
function of 70 = ^f-. The shear frequency is 52.02 Hz for (a') 
and (a"), 955.3 Hz for (b') and (b"), and 1.9689 kHz for (c') 
and (c"). 



(shadowed area in Fig. 2) is due to the invalidity of Eq. 
[5] for distances smaller than the dimension of one water 
molecule. Figure 2 indicates that the shear amplitude 
dependence of the viscoelastic modulus is very complex 
and nonlinear. For this reason we have performed mea- 
surements over a large range of shear amplitudes and fre- 
quencies (0.06 nm< X Q < 2.8 nm, 50 Hz < ui < 2 kHz). 
Following the Maxwell model for a linear viscoelastic 
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FIG. 4: The effective relaxation time as a function of 70 = 
70 • u> at d = 0.4 nm. The dash line is the fitting with Eq. [4] 



system, the relationship between the intrinsic relaxation 
time, tq, and the moduli, G' and G" , is given by 
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where Go is a constant. According to Eq. [31 G' and G" 
do not depend explicitly on 70. However, many complex 
fluids experience a drastic decrease of their structural re- 
laxation time when they are subjected to large strains. 
This phenomenon gives rise to a strong strain dependence 
of G' and G" , which can be described by the introduc- 
tion of an effective relaxation time r that depends on the 
intrinsic relaxation time and the strain rate [lfj| . Once 
defined r, it is used to replace tq in Eq. [3j and thus to 
predict G' and G" as a function of the strain. Recently, 
a phenomenological expression has been found to char- 
acterize a strain dependent effective relaxation time 



- ~ - + K ■ 70" 
r t 
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where v is a positive exponent, and if is a constant. In 
a glassy system which shows slow dynamics ( oj ^ — 



21] 



By replacing To in Eq. [3] 
— , the maximum of G" is 
Figure 3 



v ~ 1 and K ~ 1 [2 
with r in Eq. 0] when to ^> 
near 70 ~ 1, independently of the frequency, 
presents G' and G" vs. 70 for nano-confined water, ob- 
tained by applying Eq. [2] to the measured Fl and 8 at 
three different frequencies for d — 0.4 nm. G' and G" in 
Fig. 3 show remarkable behavior: (i) the peak position 
of G" is around 70 — 1 over a wide range of frequencies; 
(ii) for 70 < 1, the viscoelasticity is dominantly clastic, 
i.e., G' > G"; and (iii) G' and G" decay to zero for large 
values of 70. These features of our nano-confincd water 
system are ubiquitous in metastable complex fluids (IgT | 
and they are captured by the argument of the strain rate 
dependent effective relaxation time. Indeed, by using 
Eq. [3J and 2] the shape of the curves presented in Fig. 3 
can be fully described. 

From Eq. [31 the effective relaxation time can be pre- 
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dieted by: 



G' 1 
G" ' ~ 



(5) 



By using Eq.[5]and the experimental values of G' and G", 
r as a function of 70 for d = 0.4 nm is determined and 
shown in Fig. 4. The effective relaxation time of nano- 
confined water decreases from 40 ms to 0.7 ms when the 
strain rate increases from 14 s _1 to 6000 s -1 . The nonlin- 
earity of the relaxation time sets in when the experimen- 
tal time scale (70) is faster than the intrinsic relaxation 
time (tq). In this case, the time response can only be 
measured effectively as a function of the experimental 
time scale. 

By fitting the data in Fig. 4 with Eq. [4] we found that 
the intrinsic relaxation time was To — 0.06 ± 0.03 s for 
K = 0.95 ± 1.49 and v = 0.84 ± 0.29. These values of v 
and K are similar to those found in complex metastable 
fluids with slow dynamics 2(J 21]. The striking result is 
that the observed r and To are orders of magnitude slower 
than the relaxation time of bulk water at room temper- 
ature, which is of the order of 10 -12 s. The fact that 
confinement can drastically slow down the dynamics of 
a fluid has been previously observed in diverse systems, 
such as colloidal suspensions [22| and polymers [B| , where 
for strong fluid- wall interactions, the glass transition tem- 
perature is shifted towards high temperatures upon con- 
finement 0] • An alternative way to view this behavior is 
to consider that the confinement defines an effective tem- 
perature of the system which is lower than the canonical 
temperature. According to a previous study [13], the di- 
electric relaxation time of supercooled water confined in 
clays at 175 K is about 0.06 s similar to the relaxation 
time found in our experiments on nano-confined water at 
room temperature. Moreover, the value of the viscosity 
measured in our investigations is comparable with that of 
supercooled water at 140 K in a 100 /im radius tube [23j ]. 
A recent study has shown that the dielectric relaxation 
time of supercooled water is very sensitive to the con- 
finement 24[. For confinement lengths of the order of 1 
nm, it was found that, over a wide range of temperatures, 
the dielectric relaxation times are always longer than in 
bulk water. In our experiments, we also observe that r 
decreases with reduced confinement, i.e., with increasing 
d. Unfortunately, for d > 1 nm the lateral force becomes 
so small that we cannot measure its value precisely due 
to low signal-to-noise ratio. The only information that 
we can extract is that the intrinsic relaxation time for 
d > 1 nm is shorter than 10" 4 s. 

In conclusion, we have studied the room temperature 
viscoelastic properties of nano-confined water, finding 
a slow dynamical behavior similar to that observed in 
metastable complex fluids. By measuring the viscoelas- 
tic modulus at different frequencies and strains, we have 
found that the intrinsic relaxation time of nano-confined 
water is about 0.06 s. This value is comparable with the 



dielectric relaxation time measured in supercooled water 
at 175 K. 
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